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У робот! запропоновано новий метод реал1зацй операцй багатосливного додавання у паралельшй модели! 
обчислення. Запропонований метод базуеться на векторних операщях, що значно зменшуе юльюсть задяних 
процесор. У робот! наведено алгоритм у паралельй модел!. Анал!з складност! алгоритму показав, що у 
паралельнй модел! обчислення юмльксть однослвних операшй 4к +59 пропорщйно залежить в1д юлькост! 
задляних процесор1в к за умови, що кожен процесор виконуе однакову кльксть векторних операшй довжини 
16. Алгоритм реалзовано мовою програмування Ореп(Т, (у.1.2) та протестовано. 
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Вступление 

Микропроцессорная техника развива- 
ется очень быстро. Увеличение быстро- 
действия за счет увеличения тактовой час- 
тоты процессора или числа процессоров 
ведет к значительному увеличению потре- 
бляемой мощности устройств, что ведет к 
дорогим последующим расходам на обслу- 
живание таких устройств, управление 
процессорами, передачу данных между 
ними, написания программного обеспе- 
чения, которые сложно адаптируется при 
изменении числа задействованных процес- 
соров. Объем электроэнергии, необхо- 
димый, например, на решение задачи учета 
ошибки округления при вычислении опе- 
раций с трансвычислительной сложностью 
становится одним из факторов, который 
влияет как на выбор модели устройства 
(последовательная или параллельная), 
числа задействованных процессоров (в 
случае параллельной модели), так и 
алгоритма реализации вычислений в этой 
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модели [1-5]. 

В данной работе рассматривается 
метод, который позволяет распараллелить 
многоразрядную операцию сложения. 
Предложенный метод базируется на век- 
торных операциях, что значительно умень- 
шает число задействованных параллель- 
ных процессоров. Так, например, при па- 
раллельной обработке № элементов про- 
цессорами, которые поддерживают вектор- 
ные операции длины 16, достаточно задейст- 
вовать не более | №/16 | процессоров. 


Существует метод [1, 2], который 
позволяет сложить и -битные целые числа 
приблизительно за О(К) = =2105,К шагов, 
если задействовать К параллельных 
вычислительных узлов (К=п) и при 
условии, что одна 1-битная операция 
сложения занимает один такт (один шаг) 
процессора. Метод называется «методом с 
предсказыванием знака переноса» (сагту- 
|оокаВеа4 а@4!юоп) и этот метод включен в 
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схемотехнику почти всех современных 
компьютеров. 

При сложении двух тп -битных 
чисел на схеме потребуется О(тп)= 
=02,т+105.п шага. Далее показано, что 
данную операцию можно убыстрить, что 
актуально. 

Постановка задачи 

Даны ти -битные целые неотрица- 
тельные числа Х и У. Необходимо 
построить быстрый алгоритм вычисления 
суммы чисел Х и У в параллельной 
модели вычисления с использованием 
векторных операций. 

Пусть Х и У представляют тп- 
битные целые неотрицательные числа, 
п=кт. И пусть 2 будет их тп+1-битная 


п-1 
сумма. Число Х = УХ ‚.2") представим 


1=0 


т-1 
в виде а где Х, = о ИЕ РТ 


р=0 


1=0п-1 является т-битным словом. 


ие 1=0,п-1. тп -битные 
числа У и 1 представляются таким же 
образом. 

Далее, для удобства изложения, ти - 
битные числа Х, У, 1 и другие будем 
называть п-словными числами, каждое 
слово которого является т -битным. 

Найдем сначала сумму п-словных 
чисел Х и У без учета переноса знака 
между словами. Получим многословное 


число М = {м М т 3 где 


М, = (Х, +, т ‚ 1=0,п-1. На основании 


п-12***? 


числа М длины п слов вычислим число 
Т длины и бит, где в каждый бит 
записывается 1, , где 
ТИМ, =М"-, 
а т 
р 
Для получения числа Г сгруппируем 


биты ЕЁ, 1=0и-1, п т бит 


1 


Т= а 1 где 
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т-1 
ТЕ 2090) 
Ь=0 


1 ИМ, =" -1 
1 = им ‚1= рт-Ь, 
0 ИМ, <2"-1 


РЕКЕ 


т 


| (1) 


п 
Число Т будет состоять из К=— 
т 


слов. Использование числа Т позволяет 
спрогнозировать возникновение знака 
переноса в целой группе из т слов, на 
которые разбиваются п-словные чисел Х 
и У при многословном сложении. 

Анализ генерирования знака пере- 
носа в группах из т слов 

Рассмотрим следующее утверждение, 
которое позволяет определять знак пере- 
носа при сложении слов длины т бит без 
вычисления знака переноса, для нахож- 
дения которого обычным способом необ- 
ходимо использование дополнительного 
слова для сохранения знака переноса. 

Утверждение 1. Если числа Х,, У, 


состоящие из т бит, такие, что 
2" <Х,+У, то (М,<Х)у(М, <У), где 
М; = (Х, +7). 

Доказательство 

Далее выражение, при котором 
2" <Х,+У, будем называть условием, 


при котором возникает знак переноса. 
Необходимо рассмотреть следующие 
случаи: 


1 а - 
2 
п т ы 
2 ао или 
А аа" 


Случаи | и 3 являются простыми. В 
1-м случае знак переноса никогда не 
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риа о 508 и 
2 2 


появляется 


условие не выполняется: 
(М,<Х)\(М, <У), (2) 
где М, = (Х,+У,) =, 
Выполняется условие, при котором 


знак переноса не появляется: 
(Х,<М)^(,<М,), 6) 


где М, = (Х,+У,) =, 


В 3-м случае числа Х,, У, можно 


1 


представить следующим образом: 


Жен + ан” 
2 2 


р 2" 
где НХ, <—, НУ <—. 
ПХ ее 
В этом случае при сложении чисел 
Х,+У, всегда возникает знак переноса 
т т 


НХ, ++ Н +5->2"-1. Сумма чисел 


по модулю (Х,+У,) „= 


=(НХ, + НУ, +2") „=(НХ,+НУ)„. Так 


как НХ, < >» НУ, < г по определению, 


то после добавления НУ или НХ, к 


правым и левым частям получаем, что 


т 


НХ. НУ <5-+ НЕ =, 


т 


НУ + НХ, < +нХ, =Х.. 


Если М‚,=НХ,-+НУ, то условие 
(М, <Х,) ^(М, <7У) выполняется строже 
с условием «л», вместо «У ». 

Рассмотрим 2-й случай, когда 

5) 2" 2! 
Х, <—, —<У<2” или У<-, 

2 2 2 
о" 
т Х,<2”. Значения в числах Х,, У 
можно поменять местами, поэтому для 
упрощения можно считать, что достаточно 
рассмотреть любой из двух вариантов. 
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Выберем вариант Х, < е : = о, 
9" 
Представим число У, в виде У, = НУ + 2 
2" 
где НУ < т . Тогда 


(Х,+У), (кин 2” 
т 


Рассмотрим два случая: 


т 


а) если х.+нН <>, то знак 


переноса не появляется при 


т 


м,=х,+НЕ +5. Условие (2) не 


выполняется. Выполняется условие (3). 


т 


6) если 5х, + НУ, то знак 


переноса появляется при 


т 


м‚=х.+нЕ-5- и условие (2) 


выполняется, так как М,<У или 


т т 


Х,+Ну-5-<НУ +5, или аа 


Утверждение | доказано для всех 
случаев, что (2) выполняется, когда 
присутствует знак переноса, и не 
выполняется для остальных случаев. 

Утверждение 2. Если числа Х,, У, 
состоящие из т бит, такие, что 
(М, <ХУ(М, <У),М, = (Х,+У,) т, то 
ХИ: 

Используя Утверждение 1, на 
основании п-словного числа М получим 
п-словное число С, состоящее из вхо- 
дящих знаков переноса, которые возни- 
кают при сложении слов Х,‚ +7, 

| у (М, <ХУ(М, <У)) 
"|0 ИСХ, < М)^ < М)’ 


1=0,п-2, с =0. 
Знак 
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1 У (М, <Х,)Уу(М,,<У,) 
С = 
" 0 У(Х,,<М,,)^(,,<М,,) 


обозначим отдельно. Далее по тексту эле- 
менты, которые принимают или хранят 
знаки переноса предыдущих пословных 
слов, будем называть элементами, содер- 
жащими входящие знаки переносов, чтобы 
отличать от исходящих знаков переноса. 
Если из контекста понято, что речь идет об 
исходящем знаке переноса, то слово 
«исходящий» может опускаться. 

Для получения К-словного числа С 


сгруппируем входящие знаки переноса с,, 


1=0п-1, по т бит и запишем их в слова 
длиной т бит С ={С, _,...С)}, где 


т-— 


С,= 2:2), 
=]. {1 0<р 
Е И р=0 


Е ах и 
6. 
[0 ИХ, +У, <2" 


— 
= 
ыы 
= 
ыы 
ыы 
я 
я 


"О о ое ее 3 ВК 


РЕ рт+Ь, р= ОК, = 

т 
(С =0: (4) 

Числа С=С, |..С и ТЕТ... 


имеют одинаковую длину в ПИ бит или 


7 


Г) 
К = — слов. 
т 


Лемма 1. Если при сложении р-х 


групп ИЗ т -словных чисел 
СВ + СЕ 


однословные числа С, и Т, ‚, полученные 


по формулам (1), (4) на основе чисел 
Хх и: У}, такие, 


ртчт-1>***? ртчт-12***? 7 рт 
что С,+т. 22”, тю в ри. труше 
появляется знак переноса. 

Доказательство. 

Не уменьшая общности, доказатель- 
ство приведем для т=16. Далее на рису- 
нках х обозначает, что значение бита не 
важно. На рис. | представлен простой 
случай, когда генерируется знак переноса. 


Хх хх п х хх Хх хх 


Рис. 1. Генерация знака переноса в р-й группе из т слов. Случай 1. 


Старший бит т-битного числа С, 
равный 1 говорит о том, что сумма слов, 
предшествующих старшим словам в р-й 


У > й) т 


рт+т-2 — 


(при т=16, Х ива + Увы 22" ). т- 


группе из т слов, Хи. + 


битное число С „ содержит биты входящих 


знаков переносов поэлементных сумм т- 


словных чисел 
О }+ СА }. Слова 
Хх ИН У а Относятся к предыдущей 


группе р-1 из т слов. 


Старший бит т-битного числа Т, 
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равный 1 говорит о том, что сумма 
старших слов в р-й группе из т слов 


Х +У 


„т  Равна 2"-1 (при 
т=16,Х ив + Уже =2" -0). При до- 


рт-чт-1 


бавлении любого ненулевого значения к 
2” —1 возникает знак переноса в 
следующей группе р+1 из т слов. 

Далее на рис. 2 рассмотрим другой 
случай. 
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ооо ооооооооооогс 


т о С: ООО АВ РОО ООО 


1 


У ПО ПО Ооо ООО ПОЮ ООО № 


Р 


|2. жж алла 


х 


Хх хх п х хх Хх Хх 


Рис. 2. Генерация знака переноса в р-й группе из т слов. Случай 2 


Младший бит числа С, равный 1 
говорит о том, что при сложении слов в 
предыдущей р-1! группе из т слов 
Х +И 
знак переноса, который передается в р-ю 


группу из т слов. 
Биты числа Т, равные 1 говорят о 


(ринг + Ир, ГЕ=О,т-1, появляется 


том, что все пословные суммы в р-Й 
2" -1, 


к=0,т-1. 


группе из т слов 
Хх У =" -1, 


ртч+г ртч+г 


равны 


Очевидно, что при добавлении любого 
ненулевого значения к сумме слов 


Хт+Уиш=2" -1 приводит к тому, что 


возникает знак переноса для всей группы 
р из т =16 слов, что отражает рис. 2. 


Таким образом, если у однословных 
чисел С, и Т,, полученных на основе т - 


словных — чисел Хх а го и 
{у И — из р-й группы, согласно 


(1) и (4), какие-либо соответствующие 
биты равны 1 и все последующие старшие 
биты числа Т, также равняются 1, то это 


будет говорить о том, что в р-й группе 
при сложении генерируются знак пере- 
носа. В этом случае С,+Т, будет больше 


или равно 2”. 
Для полноты доказательства далее на 
рис. 3 рассмотрим случай, когда соответст- 


вующие биты С, и Т, не могут быть 


одновременно равными 1. 


ооо ооооооооОооОтго с 
То. фл тот 
И ЕЕ ео 


Рис. 3. Знак переноса в р -ой группе из т слов не появляется. 


Биты т-битного числа С, равные 1 
говорят о том, что при пословном 
сложении первых двух слов в группе р из 
О 
генерируются исходящие знаки переноса 
пословно, которые заносятся в 1-Й и 2-Й 
биты т-битного числа С,. 


т слов 


рт- 


Биты т-битного числа Т, равные 1 
говорят о том, что в группе р из т слов 
пословные суммы равны 
2" -1,Х =2" -1,г=3,т-1. 2-й 
бит числа Т, равный 0 говорит о том, что 


т 
Х ино = <2" -1. 


+У 


рт-г рт-г 
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В случае, когда у=1, выражение 
С.Г, 


невозможно, так как 2-й бит числа Г, 


больше или равно 2”, что 


нулевой (т.е. Х „2 +У о <2" -П). Даже, 
ай 


„а больше или 


если сумма слов Х„ 


равна 2” 
который 


и появляется знак переноса, 
добавляется к последующей 


< т 
пословной сумме Х„„, +У„,+1<2”, то 


рт+ 
этого недостаточно, чтобы появился знак 
переноса во 2-й пословной сумме. 
Покажем, что на рис. 3 значение у 
всегда будет равно 0, так как 1-й и 2-Й 
биты С, равны 1. Так как 1-й бит в С, 
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+У 


„во больше или 


равен 1, то сумма Хо 


равна 2”. Напомним, что каждый бит т- 
битного числа С, содержит входящие 
знаки переноса. Для того, чтобы нулевой 
бит в Т, равнялся | необходимо, чтобы 
выполнялось 
(Х ино + 


рт-+0 


условие, что 


И: т _ 
ВЫ я =2” -1 что невозможно 


при начальном условии, что сумма СЛОВ 
Х +У 


"о больше или равна 2” (1-Й 


рт-+0 
битв С, равен 1 на рис. 3). Так сумма по 
модулю наибольших значений равна 
О 
меньше 2” —1. 

С учетом того, что у=0, то $ =1 на 
рис. 3. Значение { зависит от суммы 


чисел Х ны Ф У р 
. т 
И 
. т г 
0 Е = 
Лемма доказана. 
Лемма | является обобщением 


ооо0оо0о0оо0о0о00 
р о Арт 


ххх хх 


функции #(.) (вепегие — генерации знака 
переноса [1]) на векторные операции. 
Лемма 2. Если при сложении р-х 


групп ИЗ т -словных чисел 
а 27° + т 2 о } 


однословные числа С, и Т, , полученные 


по формулам (1), (4) на основе чисел 
|х ии У„|, такие, 


рт-+т-12?***? рт-+т-12***?7 рт 


что Т, =2" —1, тов р -й группе возникает 
знак переноса, если знак переноса 
возникает в предыдущей р-1 группе из 
т слов. 

Доказательство. 

Не уменьшая общности, доказа- 
тельство приведем для т=16. Далее на 
рис. 4 представлен простой случай, когда 
знак переноса, полученный в р-1 группе 
из т слов передается "транзитно" в р+1 
группу из т слов. 


0 
1 


==. 


0 
1 


= 


Хх х Хх Хх Хх Хх Хх 


Рис. 4. Перенос знака из р-1 группы в р +1 группу из т слов. 


Нулевые значения бит т-битного 


числа С 


‚ говорят о том, что при 


т 
пословном сложении Х А 


рт-+г-1 за 
т=0,т-—1, т.е. знаки переноса не 
появляются. 


Биты числа Т, равные | говорят о 
том, что суммы слов равны 2”-1, 
Х Рае. Т=От- Это 
случай, когда сумма чисел 
СН У„} дает 


рт+т-12***2 рт 
число, все биты которого равны 1 я При 
добавлении знака переноса из р- 1 


группы из т слов получаем число, все 
биты которого равны 0, и при этом 
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возникает знак переноса, который 
передается в р +1 группу из т слов. 
Лемма доказана. 
Лемма 2 является обобщением 
функции р(.) (рторавзие — размножение 
знака переноса [1]) на 
операции. 

Лемма 3. Если при сложении р-х 


групп ИЗ т -словных чисел 
О 92° У. + И 29.2 У } 


однословные числа С, и Т, ‚, полученные 


по формулам (1), (4) на основе чисел 
|х ИО й У„|, такие, 


рт+т-12***2 рт 


векторные 


рт-ат-12?***? 
что С,+Т, <2", тов р-йЙ группе знак 


переноса не появляется. 
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Доказательство переноса, полученный в р-—1 группе из т 


Не уменьшая общности, доказа- слов, передается в текущую группу р из 
тельство рассмотрим для т =16. Далее на 


рис. 5 представлен случай, когда знак 


т слов. 


оооосооо оо оо ооо оотс 
ок ов т. ый Ч о 9 
0-1. Ел 
Рис. 5. Знак переноса в р группе из т слов не появляется 
Младший бит в т-битном числе С, бит числа Т, равный 0 говорит о том, что 
равный 1 говорит о том, что при сложении Х то +Уто #2" —1. Эта сумма 
слов Хил Ирлииа»  ГЕОт-И, «поглощает» знак переноса, который 
р-1! группы из т слов генерируется знак появляется в группе р-1 из т слов. 
переноса, который передается в р группу На следующем рис. 7 показано, что в 
из т слов. случае, когда соответствующие биты чисел 


Биты т-битного числа Т, равные 1 С, и Т, не равны 1, то знаки переноса 


говорят о том, что в группе р из т слов ОНА СЯ»: 


На рис.7 С,+Т, не превышает 2”. 


_ эт 
пословные суммы Х,-+Т м, =2” -1, 
Е . Лемма доказана. 
к=1,т-1. 0-й 
|1. 30 ОЕ в №0009 0 9.20: 0.20 36.56, 
а. ЗОО о Е а ЕО. 0, 
О т Ел 
Рис. 7. Знаки переноса «поглощаются» в группе р из т слов. 
Лемма 3 является обобщением знак переноса сама по себе. 
функции К(.) (Ш — «поглощение» знака На основании доказанных Лемм 
переноса ММ на векторные операции. предлагается алгоритм «предсказания» 
Примечание для Лемм 1, 2, 3. Так знаков переноса пословного сложения п- 
как в р -Й группе из т слов сумма словных чисел, с ДЛИНОЙ слова в т бит, В 
последовательной модели вычислений. 
старших слов Х„„.+Уми„. Может 


определять знак переноса в р-+1 группу из 
т слов, то считаем, что сумма 


Хх У „1 <2”-1 не генерирует 


рт+т-—1 к рт+т-— 

Алгоритм 1. «Предсказание» знака переноса при сложении чисел ДЛИНОЙ пт бит. 
п 

Вход: Числа Х ={Х,_,....Х}, У={_,...У, К= = 


п- 


Выход: С= {С,_,...С, Т= Я... М ={М,»...Мо}, с-знак переноса. 
Шаг 1. с<_0. 
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Шаг 2. Если Хи +У > 2”, то с<1. 


Шаг 3. Вычислить число М = М»... „Мо, где М; = (Хх, +, т ‚1=0,"-—1. 


Шаг 4. Вычислить число т: где Т=»(.2?), #= 


РЕ рт-Ь, РЕКЕ КЕ. 
т 


т- 1 ИМ, =2” _1 
© 0 ИМ, <2"-Г 


т-— 


Шаг 5. Вычислить число С = ото где С, = РС о 


Г (М: <Х,)У(М,, < У, ) 


га 0 М0<р 
ПТ Ир= 


с. = РЕ рт+Ь, РЕК КЕ. 


0 УМ )^а5М,)’ 
Шаг 6. Для р< 0 до р<к-2 
Шаг7. Если С,-+Т, >2", то С 
Шаг 8. Конец цикла по р. 

Шаг 9. Если С, +Т,,>2”, то с<1. 


р 


Теорема 1. Число однословных 
операций в Алгоритме 1 имеет вид 
О(п) =бп+3Е-2, где п = т -— число слов 
в каждом слагаемом, т- длина слова в 
битах. 

Доказательство. Будем считать, что 
все однословные операции сложения и 
сравнения выполняются за одинаковое 
время. Тогда Шаг 2 занимает 1 операцию. 
Для вычисления числа М на Шаге 3 
необходимо выполнить и операций. Для 
вычисления числа Т на Шаге 4 необ- 
ходимо выполнить и сравнений и столько 
же битовых операций. Если считать, что 
выполнение каждой битовой операции 
занимает столько же времени, как и 
операция однословного сравнения, то на 
выполнение Шага 4 необходимо 2п 
операции. Для вычисления числа С на 
Шаге 5 необходимо выполнить 2и 
операции однословного сравнения и п 
побитовых операций, что составляет Зи 
операции. На Шаге 7 сложение С, -Т, 
является двухсловной операцией, так как 
возможно переполнение результата сложе- 


НИЯ. Для того, чтобы уменьшить слож- 
НОСТЬ вычислений, операцию вида 
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У: 


т 


А+В>2” можно заменить следующими 
М < (А-В), 


(М <А)у(М <В), согласно доказанному 


операциями: 


ранее утверждению, что потребует 3 
операции. Таким образом, для выполнения 
цикла на Шаге 7 понадобиться З(К-2) 
операции. Шаг 9 также потребует 3 
однословных операций, аналогично 
Шагу 7. 

Общее число однословных операций, 
необходимых для выполнения Алгоритма 
| составляет 
1+1+2"+3п+3(К-2)+3=6и+3К-2. 

Теорема доказана. 

Вычисление знака переноса п- 
битных чисел С+Т дает значение знака 
переноса суммы пт -битных чисел Х и7У 
(простой случай, когда при сумме старших 
слов чисел Х и ТУ появляется знак 
переноса, не берем во внимание, так как с 
этого случая начинается выполнение 
Алгоритма 1). По сути, задача вычисления 
знака переноса при сложении ит -битных 
чисел сведена к задаче вычисления знака 
переноса суммы п-битных чисел, что в т 
раз меньше начальной длины. В свою 
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очередь, вычисления знака переноса 
суммы п-битных чисел С+Т может быть 
сведено к задаче вычисления знака 


в Г) 
переноса суммы чисел длинои в — бит и 
т 


так далее за счет рекурсивного вызова 
Алгоритма 1. Таким образом, число 
уровней рекурсивного вызова Алгоритма 1 
определяется выражением 105, п, где т — 
длина слова в битах и п -— длина числа в 
словах. 

Так, например, если п = Кт, то число 
уровней вызова Алгоритма 1 будет равно 
02, п=ю2,Кт= =1+Ю02,К. Если К<т, 
то [1+1ю2„К |=1, то рекурсивный вызов 
Алгоритма 1 не требуется. 
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Если п= Ат? и кК<т, то 
|Позп. — [Поз „би? |= [2 + © 2.К| =2, что 
говорит о том, что нужен один 
дополнительный уровень рекурсивного 
вызова Алгоритм 1. 

Вычисление знака 
пословных сумм 

Для вычисления общей суммы ит- 
битных чисел необходимо также найти 
знаки переноса пословных сумм. 


Пословные операции над С, И Г, 


позволяют вычислить знак переноса для 
следующей р-+1 группы из т слов. Для 
этого вернемся к рис. 2. и рассмотрим его 
в следующем виде. 


переноса 


0000000000000 00 1 С. 
Тов». Фе «ТЫ „О ОС В 3 ЛЬ: Ле. 1 г 
т Е Зе о ДЕ р Е 36 


Рис. 8. Корректировка пословных входящих знаков переноса С’, (в знаменателе) на основе 


входящих знаков переноса С’, (в числителе) в группе р из т=16 слов 


Младший бит в С, равный 1 говорит 
о том, что при сумме р-1 группы из т 
И 


(р-Отч+г ? 
появляется знак переноса, который 
передается в р группу из т слов и 
должен быть учтен при 
ПЕР Е 


рт+0 


слов Хх к=0,т-1, 


(р-Ют-к 


сложении 


рт+0 
Биты т-битного числа Т, равные 1 


говорят о том, что в группе р из т слов 


пословные суммы равны 
А Ин -1. г=0,т-—1. 
Очевидно, что При добавлении знака 
переноса к сумме ХУ +1 


приводит к тому, что возникает знак 
переноса во всех словах группы р из т 


слов, что показано на рис. 8. Напомним, 


что С „› р= 0, п—1, содержит входящие 
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знаки переносов. 

Для нахождения знаков переносов 
для каждого слова необходимо исполь- 
зовать побитовые операции. Условие, при 
котором соответствующие биты двух т- 
битных слов С, и Т, равны 1, говорит о 


том, что знак переноса добавляется к 


значению 2”-| и новый знак должен 
переноситься в следующее — слово. 
Предлагается следующая побитовая 
операция для итерационной корректировки 
входящих знаков переноса 


С, = С, у((С,^Т,) << 1) 


побитовая операция «<<1» сдвигает 
значения битов влево (в сторону старших 
битов) на один бит. Далее на рис. 9 
показан результат одной итерации с 
использованием побитовой операции. 


=’ где 
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ооооооооооооооотг с 
По. 3 ЗЧ ны ОЕ Ла ое ово а  РДе Е 
ооо оо бобо ооо 


Рис. 9. Корректировка входящих знаков переноса С, (в знаменателе) на основе 


входящих знаков переноса С’, (в числителе) в группе р из т=16 слов 


Для нахождения остальных бит числа 
С, необходимо выполнить дополнительно 
14 итераций. 

Вычисление пословных сумм 

Для нахождения пословных сумм с 


учетом  пословных знаков переноса 
предлагается Алгоритм 2, который 
оперирует числами С = {Сас 
ТЕТ ТИ, ММ ,..М), с, 
полученными в результате работы 
Алгоритма 1. 


Алгоритм 2. Вычисление пословных сумм 
с учетом пословных знаков переноса. 
Вход. С= Сс ь ТЕРЬ 
М={М, .,... 
Выход: 7 = Е. 

Шаг 1. Для р 0 до р<К-1 
Шаг2. Для г< 0 ог<т-1 
Шаг3. С,<С,У((С,ЛТ,) << 1) 


„Мот, с — знак переноса. 


— 
Шаг 4. Конец цикла по г. 

Шаг 5. Конец цикла по р. 

Шаг 6. Для р 0 до р<К-1 

Шаг 7. с < С,. 

Шаг 8. Для г< 0 дог<т-1 


Шаг 9. И „., < (м, + о 


рт-г 


Шаг 10. сс >>1. 


Шаг 11. Конец цикла по г. 
Шаг 12. Конец цикла по р. 

Теорема 2. Число однословных 
операций в Алгоритме 2 имеет вид 
О, (п) = 6Кт + К, где п = кт — число слов в 
каждом слагаемом, т- длина слова в 
битах. 

Доказательство. Будем считать, что 
все однословные операции сложения и 
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сравнения выполняются за одинаковое 
время. Шаг 3 занимает 3 побитовые 
операции в циклах, которые вычисляют 
знаки переноса для каждого слова. Для 
завершения циклов на Шагах 1, 2, 3, 4, 5 
необходимо Зкт побитовых операций. На 
Шаге 7 в цикле выполнится К операций. 
На Шагах 9, 10 необходимо выполнить 2 
побитовые операции и одну операцию 
сложения в циклах, которые находят 
суммы слов с учетом знаков переноса для 
каждого слова. Для завершения циклов на 
Шагах 6, 8, 9, 10, 11, 12 необходимо 
выполнить ЗАт операции. 

Общее число однословных операций, 
необходимых для выполнения Алгоритма 
2, составляет ЗАт + К + ЗАт = 6т-+ К. 

Теорема доказана. 

Циклы на Шагах 1-5 и Шагах 6-12 
могут быть объединены, так как они 
содержат одинаковые параметры для р и 
г. Циклы разделены, чтобы показать, что 
Алгоритм 2 состоит из двух логических 
частей. Первая часть вычисляет знаки 
переносов для каждого слова и вторая 
часть находит окончательные значения 
слов результата с учетом рассчитанных 
пословных знаков переноса. 

Сложение многословных чисел в 
параллельной модели вычислений. 

Разобьем п-словные числа Х и У на 
группы по 16 слов, тогда числа Х и У 
можно представить следующим образом: 
ХЕ 6 ов, бью 


ге Х16,, И!6,, р=Ок-, К=. 


16 
Далее В алгоритме используются 
векторные операции вида 


М16 < (Х16, +716) „ ‚ что равнозначно 
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выполнению цикла т; < (х, + У,) эт › 


[= 0,15 ‚ где т, , х,, у, — соответствующие 
элементы векторов М16б, Х16, : У16, ь 


Результатом операции 
М16 < (Х16, +716, будет вектор 


элементов, который будет содержать 
поэлементную сумму по модулю 2” 
(т=16) векторов Х16,, У16,. 

С16[0] - нулевой элемент в векторе 
С1б. 

С16 < С16 <<1 -— сдвиг элементов в 
векторе на одну позицию в сторону 


старших элементов: 
С160,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15] < 


С160,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14] 
016 #С!6=1 
С16 
0 7С16=0 
операцию можно переписать в виде цикла: 
М ИСИ=1 
С16[Й < й И Е ь 
0 И СП =0 
С8 < Н(С16) + [(С16) — операции 
Н(С16) и [(С16) обозначают вектора из 8 
элементов, которые содержат 8 старших и 
8 младших элементов вектора С16. 
С16П, С16, - указывают на один и 


= эту 


1 = 0,15. 


тот же элемент вектора С16б. 

Далее представлен алгоритм сложе- 
ния, который выполняется в ядре парал- 
лельного процессора, при распараллели- 
вании многословной операции сложения. 
Алгоритм 3. Ядро вычисления операции 
сложения с использованием векторных 
операций длины 16 и операций над 
словами длиной т = 16 бит. 

Вход: Числа Х = {Х16,_,....Х16,}, 


У = {716.716}, 
где К=—, р л - номер процессора, 


С16={С.,...б @.=2', 1=0/5. 
Выход: С = ее 
7 ={716,,.... 216}. 
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Шаг 1. в. — 0. 
Шаг 2. С, < 0. 
Шаг3. М16 < (Х16,+716,) в. 
Шаг 4. 

1 (М6 <х16 )\(М16<У16 
а | о. < — о < __ 
Шаг 5. С.| < С161$]. 


Шаг 6. 
С160,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15] < 


С16[0,0,/1,2,3,4,5,6,7,8,9,10,11,12,13,14] 
Шаг7. С16[0] <—0. 
С16 #С16=1 
Шаг8. С16 < | 
0 #С6=0 
И М16=2° —1 
И М16<2°-Г 


рН 


Шаг 9. 716 < . 


Шаг 10. 
С4 < (Сб) + Н([(С16)) + 
+ [(Н(С16))+ Н(Н(С16)) ` 
Шаг 11. 
Т4 < [([(Т16)) + Н([КТИ6)) + 
+ [(Н(Т16б))+ Н(Н(Т16)) ^ 
Шаг 12. 
С < (КС4))+ Н(((С4)) + 
+ [(Н(С4))+ Н(Н(СА)) 
Шаг 13. 
Т1 < [(1(Т4))+ Н(((Т4)) + 
+ (Н(Т4))+ Н(Н(Т4)) ° 
Шаг 14. Е1 < С1. 
Шаг 15. Для 1<0 до 1<К-1 
Шаг 16. Р1< РИУ(С,).. 


ее | и" =: те 
0 ИА!+Т1<2 
1 И(ЕТТ< ЕБУ(ЕТ< ТО 
| и (215 ЕГО л1< РТО’ 
ЕТ1 < (Р1+Т1) .). 


Шаг 18 ДОЖДАТЬСЯ ОСТАЛЬНЫХ 
ПРОЦЕССОРОВ 
Шаг 19. ССС 


рН рН 1етр * 
Шаг 20. Конец цикла по 1. 


(или 
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Шаг 21. ДОЖДАТЬСЯ ОСТАЛЬНЫХ 
ПРОЦЕССОРОВ 
Шаг 22. С1 < СУС.,. 
Шаг 23. Для 1<—0 до 1<16 
Шаг 24. С1< С1У ((С1лТО <<1) 
Шаг 25. Конец цикла по Г. 
1 И0< (Сл 016) 
0 /(С!лС16) =0' 
Шаг 27. 716, < С16+ Мб. 

Пояснения к Алгоритму 3 

По завершению Алгоритма 3 
= (716, ,....216,} будет содержать 
Х ={Х16, |,...Х16, 
У = 71616 а значение в элементе 
С = о, будет 
содержать знак переноса, если он возник 
при сложении чисел. 

Значения Х ={Х 16, ,....Х16}, 
У ={716, ,...,У16}, 7 = 716, .,...,716}, 
С = Пе (16 = С... находятся 
в глобальной памяти. Каждый из К 
процессоров имеет доступ ко всем 
элементам: некоторые только на чтение, 
некоторые на запись, а некоторые и на 
чтение, и на запись. 

Массив элементов 


216 * 


Шаг 26. С16 <— 


сумму чисел 


С вектора 


СЕЛО! 
используется К — процессорами для 
передачи знака переноса между собой. 
Элемент С, по завершению Алгоритма 3 


содержит знак  переноса› который 
возникает при сложении чисел 
Хх ={Х16, ,....Х16,},У = 716, |,....716 Эле 
менты — вектора С = о. С 
Р=ЬК-1 
переноса. В С, записывается 1, если при 


р ы 


содержат входящие знаки 


сложении процессором р-1 группы из 16 
слов возникает знак переноса. Процессор 
за номером р=0О обрабатывает самую 
младшую нулевую группу из 16 слов, в 
которую не передается входящий знак 
переноса, поэтому С, =0 всегда. 


Каждый р-й процессор из К 
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процессоров работает со своим входным 
вектором Х16,, У16, и результат 
сложения также записывает в свой 
выходной вектор 216, ‚ рР=Оо,К-1. 

«Областью вектора 
М16, С16, Т16 является исполняемый 
процессор. Другие процессоры оперируют 
со своими векторами. 

На Шаге 1, 2 происходит инициа- 
лизация массива, который содержит знак 
переноса. Каждый р -й процессор инициа- 


видимости» 


лизирует свой 
ЕС 

На Шаге 3 вычисляются суммы 
элементов по модулю. Каждый процессор 
работает со своей группой слов Х16б,, 
716, ‚р=о,К-1. 

На Шаге 4 вычисляются исходящие 
знаки переноса. Каждый элемент вектора 
С16 содержит знак переноса для каждого 
слова. Значения знаков переносов не 
учитывают влияние сумм элементов из 
предшествующих по номеру векторов. На 
этом этапе элементы вектора С16, С16Й, 


элемент С „ массива 


1=0,15 содержат исходящие знаки пере- 
носа пословных (поэлементных) сумм 
векторов Х16,, 716,. 


На Шаге 5 исходящий знак переноса, 
полученный при сложении самых старших 
элементов векторов Х16,, У!6,, пере- 


дается следующему процессору р+1 
посредством массива С = оо кото- 


рый содержит входящие знаки переносов. 

На Шаге 6 происходит сдвиг 
элементов, которые содержат знаки 
переносов, на одну позицию в сторону 
старших элементов. После сдвига исходя- 
щие знаки переноса становятся входящими 
и теперь можно воспользоваться Леммами 
1, 2, 3 для предсказания знаков переносов 
в каждой группе из 16 слов, которые 
представлены своими векторами М16 для 
каждого процессора. 

На Шаге 7 обнуляется элемент с 
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индексом 0 в векторе С16. 
На Шаге 8 элементы С16=С....С, 


используются только в режиме чтения и 
необходимы для того, чтобы преобра- 
зовать вектор длины 16 (С16, Т16) в слово 
(СТ,ТГ) из 16 бит, каждый бит которого 
равен 1, если соответствующий элемент 
вектора имеет ненулевое значение. Так, 
например, элементы вектора 
С6= усл е 
= {0,0,0,0,0,0,0,0,0,0,0,1,1,0,1,0} после Шага 
7 преобразуются в следующие значения 
С16 = \0,0,0,0,0,0,0,0,0,0,0,16,8,0,2,0} 
016 = {С.С} также используется 


для обратного преобразования, когда 
записываются единицы в элементы 
векторов, если соответствующие биты 
слова не равны нулю. 

На Шаге 8 знаки переносов 
преобразовываются в биты согласно 
индексам элементов вектора С16б. 

На Шаге 9 аналогичное 
преобразование происходит для вектора 
Т16. 

На Шаге 10 вектор С16 из 16 
элементов разбивается на 4 группы по 4 


слова. Все группы — складываются 
поэлементно. Так, например, если вектор 
С16 содержит элементы 


С16 = 10,0,0,0,0,0,0,0,0,0,0,16,8,0,2,0} „то 
вектор С4 будет содержать элементы со 
значениями С4 = {8,0,2,16}. 


На Шаге 11| аналогично С4 
вычисляются значения для вектора Т4. 

«Упаковка пословных входящих 
знаков переноса» в битах одного слова на 
Шагах 10, 11, 12, 13 и ‹распаковка 
пословных входящих знаков переноса» из 
битов одного слова на Шаге 26 являются 
основной идеей Алгоритма 3. После 
«упаковки пословных входящих знаков 
переноса» можно оперировать всеми 
пословными знаками переноса одной 
группы из 16 слов одновременно в одной 
однословной — операции. — Однословная 
операция сложения не применима в этом 
случае. Так, например, если входящие 


134 


155М 1561-5359. Штучний 1нтелект, 2018, № 2 


знаки переносов при пословном сложении 
равны С1б = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,/1}, 
а 716 = АЛЛАЛллАял 1}, то 
откорректированные входящие знаки 
переносов будут следующими 
Пы} . 

На Шаге 12 происходит окон- 
чательная «упаковка пословных входящих 
знаков переноса» в одном слове С1 
длиной в 16 бит. 

На Шаге 13 аналогично С1 
вычисляется Т1, который необходим для 
«транзитного» перемещения знака пере- 
носа из одного слова в следующее. 

На Шаге 14 значение СТ копируется 
в Е, который используется для вычис- 
ления исходящего знака переноса группы 
р из т слов с учетом входящего знака 
переноса из предыдущей группы р-1. 

На Шагах 15, 16, 17, 18, 19, 20 
вычисляется знак переноса, который 
появляется в группе р из т слов. 


Вычисление происходит на основе слов 
ЕТ, Т1 и знака переноса переданного из 
предыдущей — группы р-1. Число 
итераций является максимальным и 
равным К, так как для процессора с 
наибольшим индексом К -1 необходимо К 
шагов, пока знак переноса из самой 
младшей группы из т слов будет передан 
последовательно в самую старшую группу. 

На Шагах 15, 16, 17, 18, 19, 20 
выполняется вычисление исходящих зна- 
ков переносов с учетом «предсказания» 
входящих знаков переносов для каждой 
группы из т слов. 

На Шаге 16 знаки переноса, которые 
появились в группе р-1 и переданы в 


группу р через элемент С, , записываются 


25 
в самый младший бит Е1, который 
содержит копию входящих пословных 
знаков переноса. 

На Шаге 17 на каждой итерации 
вычисляется бит исходящего пословного 
знака переноса (см. рис.9). 

На Шаге 18 происходит синхро- 
низация работы всех К параллельно 
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работающих процессоров по отношению к 
чтению значений элементов массива 
С= Я = который используется для 
передачи знака переноса процессору с 
номером индекса большим на единицу. 

На Шаге 19 знак переноса, 
возникший при сложении в группе р из 
т слов, которую обрабатывает р-1, 
процессор, передается следующему за 
индексом процессору через вектор 
Се 

На Шаге 21 происходит синхрони- 
зация работы всех К параллельно работа- 
ющих процессоров по отношению к записи 
в элементы массива С = ео. 


На Шаге 22 в слово С1, которое 
содержит входящие знаки переносов 
пословных слов, в самый младший бит 
добавляется знак переноса, который 
возник в предыдущей группе р-1. С, и 


содержит значение — «предсказанного» 
исходящего знака переноса группы р-1. 
Для р-й группы этот знак переноса будет 
являться входящим знаком переноса. 

На Шагах 23, 24, 25 вычисляются 
пословные входящие знаки переносов 
внутри каждой группы р из т =16 слов с 
учетом вычисленного входящего знака 
переноса для группы р. 

На Шаге 26 происходит «распаковка 
входящих знаков переноса», которые 
необходимы для окончательного вычис- 
ления пословных сумм внутри каждой 
группы р из т = 16 слов. 

На Шаге 27 пословная сумма 
корректируется вычисленными послов- 
ными входящими знаками переноса для 
получения окончательного результата. 
Исходящие знаки переноса вычислять и 
сохранять уже нет необходимости. 

Теорема 3. Число однословных и 
векторных операций в Алгоритме 3 может 
быть представлено следующими выражениями 

О (К) =4К +59, О (К)=15, 


однословных векторных 
где К — число задействованных 
параллельных процессоров. 
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Доказательство. Считаем, что все 
однословные операции сложения, сравне- 
ния, побитовые однословные операции 
выполняются за одинаковое число тактов. 
Также считаем, что векторные операции 
сложения, сравнения, векторные побито- 
вые операции также выполняются за оди- 
наковое между собой время. 

На Шагах 1,2 выполняется 2 
однословные операции. Шаг 3 занимает 
одну векторную операцию, где длина век- 
тора равна 16. На Шаге 3 операция вычис- 
ления значения по модулю 


(х 16, +716, ) в Выполняется автомати- 


чески при оперировании элементами 
вектора с типом данных в 16 бит. Для 
выполнения Шага 4 необходимо 2 вектор- 
ные операции сравнения и одна векторная 
логическая операция. Шаги 5, 7 займут по 
одной однословной операции. На Шаге 6 
выполняется одна векторная операция, 
результатом которой является пересорти- 
рованный вектор. На Шагах 8, 9 
выполняется по одной векторной операции 
сравнения. На каждом Шаге 10, 11 выпол- 
няются по 3 векторные операции сложения 
длины 4. На каждом из Шагов 12, 13 
выполняется по 3 однословные операции 
сложения. На каждой итерации цикла на 
Шагах 15-20 выполняется по 2 логические 
однословные операции, одна однословная 
операция сравнения и одна однословная 
операция по модулю 2. Всего в цикле на 
Шагах 15-20 выполнится 4К однословных 
операций. На Шаге 22 выполняется одна 
однословная побитовая операция. На Шаге 
24 выполняется 2 однословные побитовые 
операции и одна однословная побитовая 
операция сдвига. Всего в цикле на Шагах 
23-25 выполняется 3.16 =48 однословные 
побитовые операции. На Шагах 26, 27 
выполняется одна векторная операция 
сравнения, одна векторная операция 
сложения с длиной вектора 16. 

Общее число однословных операций 
выражается 2+2+3+3+4К +1+48= =4К +59. 
Общее число векторных операций 
выражается 1+3+1+2+3+3+2=15. 
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Теорема доказана. 

Анализируя Теорему 3, видим, что 
число однословных операций не зависит 
от длины слагаемых, а зависит от числа 
процессоров, задействованных в много- 
словной операции сложения в параллель- 
ной модели вычислений. Число векторных 
операций является константой и не зависит 
от длины слагаемых и числа задейство- 
ванных процессоров. 

Выводы 

В данной работе предложен новый 
метод реализации операции многослов- 
ного сложения в параллельной модели вы- 
числений. Предложенный метод базиру- 
ется на векторных операциях, что значи- 
тельно уменьшает число задействованных 
процессоров и значительно уменьшает 
объем потребляемой мощности. В работе 
приведен алгоритм, на основании которого 
вычислительное ядро СОРИ реализует опе- 
рацию сложения в параллельной модели 
вычислений на основе векторных опера- 
ций. Анализ сложности приведенного 
алгоритма показал, что в параллельной мо- 
дели вычислений число однословных опе- 
раций 4К+59 имеет линейную зависи- 
мость от числа задействованных процес- 
соров А при условии, что процессоры 
выполняют константное число векторных 
операций длины 16. Результаты тестиро- 
вания, проведенные при реализации алго- 
ритма на языке программирования ОрепСТГ, 
версии 1.2, подтвердили корректность 
теоретических результатов. 
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апа Т аге за йслепе © "ргеЧ1с(" саггу Нах ша 
этоир оР 16 \ога$. ТБе агасе дезстЬез апа 
ог аез фе Цегайуе соггесйоп оЁ фе шри 
саггу Йаз$ ш фе \ога С ма Ы\м\м5е орега- 
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901$ изше саггу Йа УШсь аге “ргефсеа” 
ш адЧюоп оЁ Фе ргеу1юи$ эгойр$ оЁ \ог4. 
Тре аогибт пар!етеп а "рге@сйоп" оЁ 
фе саггу Нахз Базе оп опе-\ога Ш орета- 
боп$ ап4 уесог орегайоп$ Юг этопрз оЁ 16 
уога$. 'ТБе а]еоги фт 1$ 4езсге4 оп Фе у’ау 
0 Бе паретешще4 аз а соге оЁа рагаПе] рго- 
сеззог ш рагаПе| то4е] о{ сотршайоп. ТБе 
а1озаг фт 1оэлсаПу соп$15$ оЁ зеуега| раз: 
Фе гап$Рег ог Фе саггу Нах {о Ще пехЕ отоир 
айег а44т» Р1еБезё \ог4$ оЁ Ще этоир; Ше 
сасайоп оЁ Фе сгцепа С апа Т, \Шсь 
шс4ез фе "расКо" оЁ саггу Йаз5 ш Ше 
у\ога С; Цегайуе сотгесйоп о{Р сагу Йаяз 
\уВтсЬ аге “ргед се” ш адаюоп оЁ Ше рге- 
у1015$ этоирз оЁ \ог4; "апрасКе" Фе сог- 
гес{е саггу Йа25 ап ад4т?» Фет ю уае$ 
ОГ Ше \ога$ оЁ Ше этоир © её Фе Впа| 
гезай. Те и5е оЁ "расКте" / "апрасКте" фе 
саггу Нае$ оЁ \ог4$ оЁ этоир ап Ше изе оЁ 
уес‘ог орегайоп$ отеаЙу гедисез Фе питбег 
ОГ ргосеззогз$ шуоуе4. ТЬе апа]у$15 оЁ 
сотр!ехцу оЁ дезспБе4 а1еогибта $Во\у$ фай 
ш рагаПе| то4е| оЁ сотршайоп Ше патфег 
ОГ опе-41еи орегайоп$ 4ерепа$ Ппеа!у оп 
питбег оЁ ргосеззог$ чзе4 ш сасшайоп 
сопз14етие Ффаё ргосеззог$ заррош уесюг 
орегайоп$ ап@ еуегу ргосеззог ехесще$ 
сопзбапе питлбег оЁ уесбюг орегайоп$ оЁ фе 
1епэ оЁ 16. А]еогибт 15 паретещеа ие 
1апоиазе ОрепСГ (у. 1.2) апа 1е\еа. Тре 
гези14$ оЁ Пе (е5ё ргоуе оЁ Феогейса! даа. 
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